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ABSTRACT 

We examine the scalings of galactic outflows with halo mass across a suite of twenty high-resolution 
cosmological zoom galaxy simulations covering halo masses from 10 9,5 — 10 12 Mq. These simula¬ 
tions self-consistently generate outflows from the available supernova energy in a manner that suc¬ 
cessfully reproduces key galaxy observables including the stellar mass-halo mass, Tully-Fisher, and 
mass-metallicity relations. We quantify the importance of ejective feedback to setting the stellar mass 
relative to the efficiency of gas accretion and star formation. Ejective feedback is increasingly im¬ 
portant as galaxy mass decreases; we find an effective mass loading factor that scales as ^ rc 2 , with 
an amplitude and shape that is invariant with redshift. These scalings are consistent with analytic 
models for energy-driven wind, based solely on the halo potential. Recycling is common: about half 
the outflow mass across all galaxy masses is later re-accreted. The recycling timescale is typically 
~ 1 Gyr, virtually independent of halo mass. Recycled material is reaccreted farther out in the disk 
and with typically ~ 2 — 3 times more angular momentum. These results elucidate and quantify how 
the baryon cycle plausibly regulates star formation and alters the angular momentum distribution of 
disk material across the halo mass range where most of cosmic star formation occurs. 


1. INTRODUCTION 
1.1. Evidence for outflows 

Galactic outflows driven by stellar feedback have 
emerged over the past decade as an integral aspect of 
how galaxies form and evolve. Outflows are commonly 
detected in s tar-fo r ming galaxies at z 0.5 (e.g. 


Pettini et al.1 2003 

; IShaplev et al.l 200jj Weiner et al. 

2009; ISteidel et al. 

2010; Genzel et al.l 201lf Rubin et al. 
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|2013[ ). and preliminary estimates of the mass loss rate 
indicate that it is of the same order as the star forma¬ 
tion r ate ('e.g. lErbll2008l : ISteidel et ai~ll20 1 Ot [Clenzel et all 
2013). The ubiquity of outflows during this cosmic pe¬ 
riod and th eir lack of a ssocia tion with active galactic nu¬ 
clei (e.g. iShaplev et al 1 [20031 ) suggests that they are not 
associated with a particular phase of galaxy formation. 
Instead, it appears that these galactic outflows are com¬ 
mon in galaxies lying on the star forming galaxy main 
sequence and are likely driven by their high star forma¬ 
tion rates or some related property. 
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Theoretical models of galaxy formation have also 
found it increasingly necessary to invoke stellar feedback- 
driven galactic outflows. Strong and ubiquitous out¬ 
flows first appeared in galaxy formation simulations in 
order to explain observations of widespread metal en¬ 
richment in__the_Jntergalacticjnedium_(IGM}at early 
epochs (1 Aguirre et al.l 120011 iQppenheimer et al.l (20061) . 
However, it was soon recognized that these same outflows 
had a considerable impact on the properties of galaxies 
themselves, such as their stellar and baryonic content, 
their metallicity, their mass distribution, and the state of 
the circum-galactic medium (CGM). Today, the majority 
of successful galaxy formatio n models, be they hyd rody- 
namic or semi-analytic (e.g. iSomerville et al.ll2015h . in¬ 
clude stellar feedback-driven galactic outflows as a cen¬ 
tral feature. 


1.2. Impact of outflows on galaxies 

Galactic outflows have a strong impact on the growth 
of the stellar and gaseous content of galaxies. First, 
mass loss limits the amount of baryons in the disk of 
galaxies, thereby impacting the observed b aryonic mass 
fraction within the halos of galaxies (e.g. Ifwl (2009) 
and the total stellar and dis k gas mass of the galaxies 
fe.g. lScannapieco et al.l (2012). Mass loss, in combination 
wit h heati ng o f the interstellar media by feedback (e.g. 
iBrooks et al.1 120071) . also impacts the observed galaxy 
stellar mass function by limiting the amount o f star for¬ 
mation in the disk of th e galaxy (|Snringel fc Hernauistl 
l2003bt IDave et aLll201l ) . By reducing the amount of gas 
available for star formation, outflows mitigate the over¬ 
cooling problem wher ein too large a fraction of baryons 
cond e nses into stars (iWhite fe Ilees]119781: iBalogh et al.l 
1200 H IDave et HI 120011) . Moreover, stellar feedback- 
driven outflows are frequently invoked to explain the 
sub-L * portion of the stellar nrass-halo relation (e.g. 
iShankar et al.ll2006l) and the simulations most successful 
in reproducing the stellar mass-halo mass relation typi¬ 
cally employ strong stellar feedback dStinson et al.ll2012l : 
lAumer et al.l l2013t iHookins et~aD 1201.31 ) . The efficiency 
of outflows across a range of galaxy masses is therefore 
key to determining the stellar and baryonic mass for dif¬ 
ferent mass halos. 

In well as globally impacting the baryons within 
galaxies, outflows play an important role in determining 
the distribution of matter. For example, the stellar 
feedback that drives outflows delays star formation, 
which results in less angular momentum loss during 
galaxy mergers a nd, therefore, le ss ce nt rally concen¬ 
trate d galaxies (I White fe Reesl 119781: iDekel fe Silk 
1986 : iNavarro fe Steinmetzl 119971 iRobertson et al. 
2004 : lO kamoto et al.1 120051: IGovernato et, al.l 120091 : 


Scannapieco et al.ll2008( ). 


In addition to being an agent for reducing angular 
momentum loss, stellar feedback is now also consid¬ 
ered a mechanism for removing low-angular momen¬ 
tum baryons in galaxies. This removal is necessary 
for the modeling of realistic galaxies: while the over¬ 
all specific angular momentum of disk s is comparab le to 
that o f their parent dark mat t er halos (iFalffc Efsta thioul 
119801 : iDalcanton et al.1 Il997t IMo et al.l I1998D . the ob- 
served distribution of the specific angular momentum 
in any given disk differs from predictions based on 
pure angular momentum conservation. In particular, 


observed galaxies are deficient in low-angular momen¬ 
tum materi a l compared to sim ple disk collapse mod¬ 
els (Ivan den Bosch et al.l 120011 ). Strong galactic winds 
were first posited as a mech anism fo r removing low- 
angular momentum material in lBinnev et all (1 2001 ) and 
recent simulations have demonstrated that stellar feed¬ 
back preferentially r emoves low-angular momentum ma- 
terial from galaxies (IGovernato et al.1 120101: iBrook et al.l 
l20T~il 1201 2H IMa.ccid et all (2101 ITThter et a,1 .1120141) . The 
loss of low-angular m omentum baryons resu lts in bul¬ 
geless dwarf galaxies (|Oovernato et al.ll2010l ). and spi¬ 
ral galaxie s with more realistic central baryonic dis¬ 
tributions (IBrook et al.1 12012at iChristensen et al.ll2014al : 
lAngles-Alcdzar et al.ll201~lb 

While outflows are typically thought of as a way to re¬ 
move material, the reaccretion of that material (wind re¬ 
cycling) also impacts galaxy evolution. Some outflowing 
material must be transferred to th e diffuse IGM in order 
to explain obser vations (e.g. iCowie et~ah I1995D of metal¬ 
line absorption (fO ppenh eimer et al.ll2006 : lCen fc Chisaril 
120111 : IQppenheimer et al.l l2012j h However, in many 
cases the outflows are thought not to escape the CGM, 
but to instead return to the galaxy on relatively short 
timescales. This so-called wind recycling adds to the 
pristine accretion from the IGM and to the accretion of 
gas already bound in galaxies (i.e. mergers). The recy¬ 
cling of previously ejected wind material can be a key 
factor in setting the galaxy stellar mass fu nction (e.g. 
IQppenheimer et aD 120101 : iBower et al.l I2012D . It is also 
thought that fountaining gas gains angular momentum 
through interactions within the halo environment be¬ 
fore being reaccreted dMarasco et al.l [2012T) . which fur¬ 
ther shifts the angular momentum distr i butio n of the 
disk baryons to higher values (IBrook et al.ll2012a ). 

1.3. The analysis of outflows presented here 

Together, inflows, outflows, and wind recycling gov¬ 
ern m any of the key physical properties of galaxies (e.g. 
IDave et aLll201 2i). Therefore, it is critical to understand 
the operation of the baryon cycle, including the scaling 
of the mass loading factor with galactic properties, the 
relative rates of outflowing, accreted, and recycled gas, 
and the source and eventual destination of outflowing 
material. Hydrodynamic simulations of galaxy growth 
are a valuable tool for this, because the inherently dy¬ 
namical nature of the baryon cycle requires a fully dy¬ 
namical model to capture it properly. Such simulations 
must be set within a fully cosmological context, since the 
accretion is cosmologically-driven. Furthermore the de¬ 
tails of the complex interactions between inflows, stellar 
feedback, outflows, and ambient halo gas strongly moti¬ 
vates very high numerical resolution. These requirements 
are a challenge for current galaxy formation models, one 
that is only recently starting to be met using cosmolog¬ 
ical “zoom” simulations, in which an individual galaxy 
is re-simulated at much higher resolution within a larger 
cosmological volume. 

In this paper we investigate the detailed dynamics of 
inflow and outflow processes using a suite of cosmologi¬ 
cal zoom simulations. We take advantage of the particle- 
based nature of our Smoothed Particle Hydrodynamics 
(SPH) simulations run with the Gasoline code to di¬ 
rectly track all mass movement in and out of the disk 
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and CGM. This particle tracking enables us to directly 
study recycling and to identify the source and future 
trajectory of individual parcels of gas. We analyze high- 
resolution simulations of twenty spiral and dwarf galaxies 
that span two and a half orders of magnitudes in virial 
mass, all simulated with the same physics and compa¬ 
rable numerical resolution. Our work improves on pre¬ 
vious particle tracking analyses that have focused either 
on lo w-resol ution (non-zoom) simulations of many galax¬ 
ies (|Onpenheimer et al.||201 0D or on a few s imulations of 
similar-mass galaxies ((Brook et al.l 1201 2at iWoods et al.1 
12014 lUbler et~ahl 120141 and com plements non-particl e 
tracking studies of galactic winds (IMuratov et al.l[20T5h . 

Our simulations self-consistently generate outflows 
from the available supernova energy using a methodol¬ 
ogy that has been shown to successfully reproduce a wide 
variety of galaxy observations. In this method, the trans¬ 
fer of stellar feedback energy depends only on the local 
properties of the ISM, and since the feedback model is ig¬ 
norant of the host galaxy properties, outflow trends with 
mass result from the dynamics of the simulation. This 
analysis enables us to independently study inflows, out¬ 
flows, and recycling as a function of galaxy mass, which 
allows a deeper investigation into the underlying physi¬ 
cal processes that govern the baryon cycle. The details of 
baryon cycling presented here do depend on our method¬ 
ology for driving outflows, and may vary f or different 
choices of physical models (see e.g. iKeller et al.l 1201, ^ 
IMuratov et al.112015} ). Nonetheless, our results are of par¬ 
ticular interest since our outflow driving model yields a 
viable match to numerous key observational constraints, 
as we will show. 

Our paper is organized as follows. In (J2]we describe the 
simulations used and the particle tracking analysis used 
to determine outflows. We justify our simulation models 
in j jU by comparing their global properties to observed 
trends. In ! ro rT~n we analyze the efficiency of various 
forms of feedback across galaxy mass. We examine gas 
recycling and the characteristics of reaccreted material 
mM m We discuss our results in light of other models 
and numerical concerns in 1|U 

2. SIMULATION AND ANALYSIS 

We compared the properties of outflows across a 
set of twenty field galaxies with final virial masses 
between 10 93 to 10 12 M 0 that were si mulated using 
the i V-Body + SPH code, gasoline (jWadslev et al’l 
120041. GASO LINE is an SPH extension to PKDGRAV 
(|Stadell l2001h . a parallel, gravity-tree based TV-Body 
Code. The simulations were integrated to a redshift 
of zero in a_ fully-cosmologica l, ACDM context using 
WMAP3 (IS per gel et al.l 120071 ) parameters: H 0 =0.24, 
A=0.76, h=0.73, cr s =0.77. 

All twenty galaxies are central galaxies selected from 
a set of eight simulations. In order to achieve signifi¬ 
cantly higher resolution while still modeling the effects 
of the large-scale environment, w e used t he “zoo m-in” 
volume renormalization technique (|Katz fe White 1993). 
More specifically, to create the initial conditions for these 
simulations we selected eight field-like regions from two 
uniform dark matter-only simulations, one representing 
a 25 3 Mpc 3 and the other a 50 3 Mpc 3 volume. These 
regions were then resimulated at higher resolution. The 
final sample of galaxies has force spline softening lengths 


of either e = 87 or 170 pc. The particle masses for 
the dark matter, gas, and stars (at their formation) 
are, respectively, 1.6 (13)xl0 4 , 3.3 (27.0) xlO 3 , and 1.0 
(8.0) x 1O 3 M 0 . The simulations have a minimum smooth¬ 
ing length of O.le, which is sufficient to resolve the disks 
of galaxies and the giant molecular clouds within which 
stars form. The parameters used to simulate the galaxies 
and their final properties at z = 0 are listed in Table 1. 
All the simulations used a force accuracy criterion of 9 = 
0.725 and a Courant condition of 77q=0.4. Particle time 


steps were required to satisfy At = p 


where r?,- is 
a t 11 


the particle’s gravitational softening, cq is the particle’s 
acceleration, and 77 = 0 .195. 

We integrate over the H and He chemical networks to 
produc e non-equilibrium ion a bundances and H 2 abun¬ 
dance ((Christensen et alJl2012l ). H 2 forms both on dust 
grains, assuming a fix ed dust-to-metallici t y ratio and a 
clumping factor of 10 (|Wolfire et al.l 120081 : iGnedin et al.1 
2009) and via H , following the minimal model of 
Abel et all (|1997l ). Photo-ionization and heating rates 
of H and He are calculated assuming a set cosmic ultra¬ 
violet (UV) background. This cosmic UV background is 
redshift-d ependent, following an up dated version of that 
laid out in lHaardt &; Madauf (11996ft . The Lyman-Werner 
radiation, which is responsible for H 2 photodissociation, 
is calc ulated based on emissio n from nearby stellar par¬ 
ticles (|Christensen et alll2012D . H 2 is shielded from dis¬ 
sociating r adiation both through se l f-shielding and dust 
shiel ding (IDraine fc Bertoldil 1 19961 IGlover fc Mac Lowl 
120071: IGnedin et al.ll2OQ90 7 using the smoothing lengths of 
particles for the column lengths. Similarly, HI is shielded 
from photo-ionizing radiation by dust. 

Cooling chan nels include collisional ionization 
( AbeLeiLaDll 997D , H 2 collisions, radia tive recombination 
( Blackl 119811: IVerner fe Ferlandl [l996h . photoioniza tion 
bremsstrahlung, and HI, H 2 and He line cooling dCenl 
1992). Addition al cooling takes place via metal lines 
(|Shen et al.ll2010l) . The metal line cooling rates used in 
the code are tabulated based on the gas temperature, 
density, and metallicity and the cosmic UV background. 
These rates are com puted using CLOUDY (version 07.02; 
iFerland et all '1 998 1 under the assumptions that the gas 
was in ionization equilibrium and optically thin to UV 
radiation. 

GASOLINE separately follows both the oxygen and iron 
abundances of gas particles and the total metal produc¬ 
tion. These metals are injected i nto th e gas by type I 
and II SN following iRaiteri et al.l (11996ft and distributed 
across the smoothing sphere. They are also injected by 


ste llar winds using a model for mass loss that follows 
IWeidemannl (1198711 and assumes that the metallicity is 
that of the stellar particle. Metal s are furthe r distr ibuted 
throughout the gas by diffusion (iShen et al.ll2010 ). 

Star formation proceedes probabilistically, according 
to the gas density and H 2 fraction. The dependency 
on H 2 represents the observed connection between the 
sta r format io n ra te and the local H 2 abundance (e.g. 
iBigiel et~ani2008l ). The probability p of a given gas par¬ 
ticle forming a star is 


p — m 9 as _ e _c * * H2 +^ H i At / t<! « n ) 
TYlstar 


(i) 
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TABLE 1 

Properties of the a set of galaxies with different ISM models at z = 0 


Simulation 

Softening 

Gas Particle 

Halo ID 

Virial Mass 

Gas Mass 

Cold Gas 

Stellar 

Vf 

Name 

Length 

Mass 



in Rvir 

Mass 

Mass 


[pc] 

[M 0 ] 


[M©] 

[M 0 ] 

[M©] 

[M©] 

km/s] 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

h799 

87 

3.3 x 10 3 

12 , 3,5 

2.4 x 10 lu 

1.4 x 10 9 

2.5 x 10 8 

1.4 x 10 8 

55 




4 

6.8 x 10 9 

4.1 x 10 7 

1.0 x 10 7 

1.8 x 10 7 

33 




6 

4.4 x 10 9 

3.9 x 10 7 

2.6 x 10 7 

3.5 x 10® 

27 

h516 

87 

3.3 x 10 3 

^ 1 , 2 , 3,0 

3.8 x 10 1U 

2.3 x 10 9 

5.5 x 10 8 

2.5 x 10 8 

67 




2 

1.5 x 10 10 

3.7 x 10 8 

4.6 x 10 7 

8.1 x 10 7 

34 

h986 

170 

2.7 x 10 4 

-^ 3,5 

1.9 x 10 11 

1.7 x 10 iu 

3.5 x 10 9 

4.5 x 10 9 

103 




2 

5.9 x 10 10 

3.2 x 10 9 

7.4 x 10 8 

1.2 x 10 9 

77 




3 

3.8 x 10 10 

2.4 x 10 9 

5.4 x 10 8 

4.6 x 10 s 

76 




8 

3.8 x 10 10 

6.4 x 10 7 

1.3 x 10 7 

4.0 x 10 7 

35 




15 

4.4 x 10 9 

8.7 x 10 7 

2.7 x 10 8 

6.2 x 10® 

29 




16 

3.2 x 10 9 

3.0 x 10 7 

1.1 x 10 8 

2.3 x 10® 

27 

h603 

170 

2.7 x 10 4 

^ 3,5 

3.4 x 10 11 

3.1 x 10 iu 

4.2 x 10 9 

7.8 x 10 9 

115 




2 3 

1.0 x 10 11 

6.1 x 10 9 

7.8 x 10 8 

3.8 x 10 9 

75 




3 

2.9 x 10 10 

1.8 x 10 8 

1.8 x 10 8 

3.9 x 10 s 

50 

h258 

170 

2.7 x 10 4 

l 3 ' 4 

7.7 x 10 11 

5.6 x 10 iu 

5.7 x 10 9 

4.5 x 10 1U 

182 




4 

1.1 x 10 10 

1.4 x 10 8 

6.2 x 10 7 

5.9 x 10 7 

43 

h285 

170 

2.7 x 10 4 

l 3 

8.8 x 10 11 

6.3 x 10 iu 

8.5 x 10 9 

4.6 x 10 1U 

164 




4 

3.4 x 10 10 

1.2 x 10 9 

1.5 x 10 8 

3.9 x 10 s 

64 




9 

1.2 x 10 10 

3.1 x 10 8 

1.3 x 10 8 

5.4 x 10 7 

52 

h239 

170 

2.7 x 10 4 

l 3 

6.8 x 10 11 

8.1 x 10 iu 

6.2 x 10 9 

4.5 x 10 iu 

165 


1 Appears in IChristensen et all J20 1 21 ). 

2 Appears in IGovernato et alJ 1201 21 ). 

3 Appears inlMunshi ct al. (20111). 

4 Appears in IZolotov et all 1120121 1. 

5 Appears in Christensen e t. all 120 1 4bl ). 

where m gas is the mass of the gas particle, m s t a r is the 
mass of the potential star, c* = 0.1 is the star forming 
efficiency, Xh 2 and Xhi are the mass fractions of the 
particle in the form of H 2 and HI, respectively, At is the 
time step, and td yn is the dynamical time. Star formation 
is allowed to proceed only in gas particles that are denser 
than 0.1 amu cm' 3 and colder than 10 3 K. However, 
the dependency on the H 2 abundance made these two 
constraints largely irrelevant since almost all stars form 
at gas densities greater than 10 amu cm -3 . 

Supernova feed back is incorp o rate d using the “blast- 
wave” approach (jStinson et alJl2006[ ), in which the the¬ 
oretical solution to a blastwave explosion in a medium 
of a given density and pressure are used to determine 
the spatial extent of the feedback. In this approach, en¬ 
ergy is distributed to nearby gas particles that h e wi thin 
t he m aximum radius of the SN blastwave (|Clievalierl 
11974 1. The cooling of these affected particles is dis¬ 
abled for a period of time equal to the theoretical life¬ 
time of the hot, low-density shell produced during the 
momentum-conserv i ng p hase of the supernova remnant 
dMcKee fc Ostrikerl I1977H . Typically, these periods of 
time are on the order of several times 10 7 years. The total 
amount of energy deposited in the ISM is the canonical 
10 51 ergs per supernova. 

The blastwave recipe di ffers from man y ot h er sub¬ 
grid feedback rec i pes (e.g. iSpringel fc Hernauistl [2003a; 
iDave et all 120111 : IScannapieco et all l2012t l in that the 
feedback-affected gas particles remain hydrodynamically 
coupled to the rest of the simulation. Additionally, the 
feedback depends only on the local gas properties; it is 
independent from the large-scale properties of the galac¬ 


tic halo. This feedback recipe does not include a sepa¬ 
rate model for other forms of stellar feedback, such as 
radiation pressure, which may help drive a galaxy wind 
either by adding additional momentum to the gas, or 
causi ng the gas to be more responsive to super nova feed¬ 
back (|Honkins et, al.ll20lil iStinson et al. 1120121) . As such, 
this efficient transfer of supernova energy into the ISM 
and the temporary delay of cooling is best interpreted as 
a model for the total stellar feedback from all processes 
related to young stars. In § we discuss the possible 
impacts of excluding other forms of stellar feedback. 

2.1. Post-processing Analysis 

Individual halos are iden tified in each snapsho t us- 
ing AMIGA’S halo finder (iKnollmann fe Knebel [20091 : 
IQill et alJl200l FI. which uses a grid hierarchy to iden¬ 
tify areas of over-denisity and iteratively removes gravi¬ 
tationally unbound particles from the prospective halos. 
The virial radius, R V i r , is defined to be the radius for 
which the average halo density is some multiple of the 
background density. This value for the average halo den¬ 
sity evolves with redshift but is approximately equal to 
100 times the critical density. In determining the evolu¬ 
tion of the halo, we use a merger tree to trace the main 
progenitor back in time. At each snapshot, the main pro¬ 
genitor is defined to be the halo in the previous step that 
contains the majority of the particles in the current halo. 

In order to verify the properties of the simulated galax¬ 
ies, we make as direct a comparison to observations as 

1 AMIGA’S HALO finder is available for download at 
http://popia.ft.uam.es/AHF/Download, html 
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possible (e.g. see S I3.1l for the comparisons to the stellar 
mass-halo mass, Tully-Fisher and mass-metallicity rela¬ 
tionships). This comparison requires the generation of 
mock-photometric magnitudes, which was accomplished 
using th e ray traci ng radiative transfer program, SUN¬ 
RISE (l.Tonssonl 12006 1. SUNRlSE-generated images of all 
the galaxies in SDSS g , r, and i filters are shown in Fig¬ 
ure u 

Table 1 lists the global properties of the halos, includ¬ 
ing virial, total gas, cold gas, and stellar mass for each 
of the halos at a redshift of zero. We list both the total 
gas mass (mass of gas particles within R V i r ) and cold gas 
mass in order to distinguish between the entirety of the 
gas mass and the fraction of it that is easily observable. 
Here the cold gas mass was defined to be gas mass in the 
form of HI, H 2 , and Hel. 

2.1.1. Particle Tracking 

In order to study the baryon cycle within galaxies in 
detail, the gas must be followed as it is accreted to and is 
ejected from the disk. As such, we determine gas accre¬ 
tion and gas loss through particle tracking. Essentially, 
we divide the gas between disk, halo, and IGM for each 
of the snapshots, and then the movement between these 
phases marks inflow and outflow. In addition to allowing 
us to identify gas recycling, this particle tracking enables 
us to determine the source of gas outflows. 

We used particle tracking to identify accreting and out¬ 
flowing gas and further subdivided the outflowing gas 
into ejected gas that became dynamically unbound from 
the disk and expelled gas that escaped the halo. Our 
method for identifying these incidents is as follows. First, 
at each snapshot we determine the gas particles that are 
in the main halo and the disk of the galaxy. Gas parti¬ 
cles are considered part of the galaxy if the halo finder 
determined them to be a member of either the main halo 
or one of its satellites. Gas particles are considered part 
of the disk of the galaxy during a snapshot if they meet 
all of the following criteria: 

• has a density n > 0.1 amu cm -3 

• has a temperature T < 1.2 x 10 4 K 

• are within 3 kpc of the midplane of the galaxy. 

The density and temperature cuts select for the cool in¬ 
terstellar medium phase of gas, while the spatial cut ex¬ 
cludes the ISM of satellite galaxies. For this analysis, we 
focus on tracking particles since 2 = 3. 

Instances of gas accretion includes both the first snap¬ 
shot a gas particle is identified as part of the disk, and 
each time it reenters the disk. 

Ejected gas is defined to be the gas particles which 
not only stop being identified as part of the disk but 
also become gravitationally unbound from the baryonic 
disk of the galaxy. These ejected particles have, at some 
snapshot subsequent to leaving the disk, a kinetic energy 
greater than the gravitational potential from the com¬ 
bined disk gas plus stellar mass (calculated as if all the 
mass were located at the center of mass of the galaxy). 
Note that these particles need not be unbound from the 
greater dark matter potential well. Gas particles may be 
ejected multiple times in their history provided they are 
identified as part of the disk in between ejection events. 


While ejected gas particles need only become unbound 
from the baryonic disk, a subset of the ejected gas even¬ 
tually also escapes the halo. Any gas particle that ends 
up beyond I? v ir at a snapshot subsequent to ejection is 
identified as having been expelled from the halo. 

In addition to the ejected gas, there is a larger popu¬ 
lation of gas particles that are labeled part of the disk 
in one snapshot, but not in a subsequent one. This oc¬ 
curs whenever a gas particle becomes too hot, too low 
density or too distant from the mid-plane of the galaxy. 
As such, it can include tidally stripped material, in ad¬ 
dition to gas affected by feedback. Ejected gas particles, 
therefore, are a subsection of this broader category just 
as the gas particles that are expelled from the halo are a 
subsection of the ejected gas. 

These identifications of gas outflow and accretion are 
limited by the temporal spacing of the snapshots. For the 
simulations in this paper, snapshots were spaced approx¬ 
imately 100 Myrs apart. Our time resolution for tracing 
inflows and outflows, therefore, is also about 100 Myr. 
Furthermore, the limited number of snapshots mean that 
particles may have been ejected and reaccreted between 
two concurrent snapshots. In such a case, an outflow 
would not have been identified. This limited time res¬ 
olution means that the outflowing and accretion rates 
must be considered lower limits. The similar snapshot 
spacing across simulations ensures, however, that we can 
draw a comparisons across the different galaxies. 

Figure [2] shows the cumulative mass loss history for 
each of our simulated galaxies versus time since the Big 
Bang, as a fraction of the final virial mass (listed in the 
upper left). The dashed line shows the total baryonic 
mass in the disk as a fraction of the redshift zero halo 
mass. Merger events can be identified by the sudden 
jumps, more common in larger-mass halos and at earlier 
epochs. 

The black line shows the cumulative gas mass accreted 
to the disk, including reaccretion events (stars are not 
included). The green line shows the cumulative outflow 
mass from the disk. While disk masses tend to stabilize 
at later epochs, the total mass loss generally increases 
roughly linearly with time, mimicking the accreted mass. 
According to these accretion and outflow definitions, a 
roughly constant mass of gas enters and leaves the disks 
at late times. Also notable is the greater amounts of 
gas accretion and loss from the disk compared to the 
final virial mass in more massive galaxies. The smaller 
baryonic fraction of dwarf galaxies is the result of both 
lower accretion rates compared to their virial mass and 
the loss of a greater fraction of their accreted material, 
shown further in § 13.21 

We further subdivide the cumulative mass loss into gas 
ejected from the disk (red) and gas that is expelled from 
the halo (blue). Since gas expelled from the halo is a 
subset of gas ejected from the disk, it is always lower, 
but the two generally track each other well. At late times 
much of the disk mass loss is not in “ejected” material 
(i.e., gas that becomes dynamically unbound from the 
disk). Rather, it is dominated by gas that escapes the 
disk but does not become dynamically unbound. This 
gas that is removed from the disk but not ejected includes 
material heated by supernovae to above 10 4 K, gas that is 
tidally stripped, and gas that is entrained by outflowing 
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Fig. 1.— Simulated observations of the sample of galaxies in SDSS r, g , and i bands at z = 0 and ranked by mass. All galaxies are 
shown face-on and the images were generated using SUNRISE. 
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3. RESULTS 

3.1. Global Galaxy Properties 

The balance between gas accretion, star formation, 
and outflows determines the baryonic content of galax¬ 
ies. Therefore, the observed stellar and cold gas masses 
of galaxies act as a basic constraint on theoretical mod¬ 
els. Our first step, therefore, is to verify our simulated 
galaxies’ agreement with global z = 0 trends in the bary¬ 
onic content of observed galaxies, in particular the stellar 
mass-halo mass, Tully-Fisher, and mass-metallicity rela¬ 
tions. 

The halos in our sample are shown to match the 
z = 0 stellar-mass halo-mass relation i nferred from 
abundance matching (iBehroozi et ail 1 201. 'll iMoster et al.1 
1201, It iBrook et al.ll2014h in Figurelfll In order to make an 
accurate comparison, stellar masses were calculated from 
simulated photometric observations while halo masses 
were t aken from d ark matter-only simulations, as in 
iMunshi et al. (j2013t ). The agreement between the sim¬ 
ulations and the abundance matching models indicates 
that the simulated galaxies are able to form the correct 
mass of stars for their halo mass over the course of their 
lifetime. 

Figure [4] compares the simulations to the ob served 
baryonic Tully-Fisher relationship from lMcGaughl (!2005l l 
(the sum of the stellar and cold gas mass of galaxies as 
a function of their rotational velocities). To ensure a 
fair comparison to these data, we compute the baryonic 
mass as the sum of the stellar mass calculated from the 
galaxies’ B-band magnitudes using a m ass-to -light ratio 
determined by the B-V color followine fMcGaugh (2005), 
together with the HI, H 2 and Hel gas mass. Addition¬ 
ally, the asymptotic velocity, Vf, is calculated by fitting 
either an increasing or decreasing arctangent function to 
the circular velocities for radii that lay between twice 
the softening and less than the radius containing 90% of 
the cold gas. The agreement is excellent over the entire 
range of overlapping masses. 

Analogously, Figure [5] shows the simulated galaxies 
along the observed SPSS A band Tully-Fisher r e lation ¬ 
ship from iGeha et al.l (120061) and iPizatnio et all (12007ft . 
The lPizaeno et alX i 2007ft data are for the tangental rota - 
tional velocities, while the data from IGeha et al.1 ([2006ft 
employs the HI line width and thus should be interpreted 
as lower limits to the circular velocities. The simulated 
galaxies follow both the baryonic and standard Tully 
Fisher relations across nearly an order of magnitude ve¬ 
locity range. 

The Tully-Fisher plots confirm that the feedback 
mechanism in these simulations is able to solve the 
so-called angular momentum crisis in galaxies noted 
in early simulatio ns wi thout str on g ste llar-driven feed- 
back (lNavarro fe Steinmetzl 120001 : iSteinmetz fc Navarrol 
HMjft . Feedback primarily solves this problem by the 
stro ng supp ressi on of star formation in dwarf galax¬ 
ies ([Munshi et alii 2l)13 i. As a result of this sup¬ 
pression, there is less growth via merging that would 
otherw ise generate an overly peaked central rotation 
curve (|Governato et al.1 12009ft . Additionally, as will 
be further discussed in these outflows reduce the 

amount of material in the centers of galaxies through the 


preferential rem oval of lotwangular momentum gas. As 
noted bv IBrook et al.1 (l2012aft . outflows play a key role 
in setting the angular momentum distribution in galax¬ 
ies, and the agreement with the baryonic and i-band T-F 
relation suggests that the outflows in these simulations 
are plausible. 

The gas-phase metallicities of galaxies are controlled 
by the balance between the accretion of pristine gas, 
the injection of metals into the ISM from stars, the 
ejection of metal-enriched ISM g as by feedback an d the 
re-accretion of enriched gas ie.g. iFinlator fc Davell2008l : 
iPeenles fc Shankar! 1201 llTPave et al.ll 2012 ft . Despite the 
complexity of the processes, ga laxies follow a well-de fined 
mass-metallicity relationship dTremonti et al.l [2004ft . al¬ 
beit subjec t to calibration uncerta inties in metallicity 
indicators dKewlev fc Ellisonl [2008). Thus the mass- 
metallicity relationship provides another strong test 
of the plausibility of our outflow model. Previously, 
iBrooks et al.l (12007ft found galaxies generated using an 
earlier version of GASOLINE follow the observed mass- 
metallicity relationship. Here, we update that analysis to 
verify our sample of simulated galaxies with the current 
ISM model and star formation and feedback parameters. 

Figure [ 6 ] shows our simulated galaxies in relation 
to the ob served redshift ze r o mass-metallici t y rela tion¬ 
ship from ILee et all (2QQp, lTremont i_et_aLj_( 20041) and 
lAndrews fc Martin i (2013)0. The lTremonti et all (12004ft 
metallicity data were lowered by 0.26 dex in order to ac- 
co unt for the o ffset i n the metall icity calibration, as noted 
bv lErb et al.1 (|2006ft . Following ILee et all (12006ft . stellar 
masses for the simulations are calculated from the K and 
B-band magnitudes. Gas particle oxygen abundances are 
weighted by the particle’s star formation rate (i.e. prob¬ 
ability of star formation) to mimic the measurement of 
metallicities in star forming regions for observed galax¬ 
ies. The simulated galaxies with stellar masses between 
10 s and 10 9 Mq have somewhat lower metallicites than 
the observed galaxies. Nevertheless, the simulated galax¬ 
ies broadly follow a power law with [0/H]oc M 9 3 from 


M* « 10 6 ’ 5 — 1 O 1O 5 M 0 , which is in good agreement with 
the observed galaxies over that wide range of masses. 

Taken together, the stellar mass-halo mass, Tully- 
Fisher and mass-metallicity relations relation provide 
stringent constraints on the cumulative effects of gas in¬ 
flows, gas outflows, and star formation. The agreement 
of the simulations with these relations imply that such 
processes and their scalings with mass are plausibly rep¬ 
resented in our simulations. We note that there has been 
no specific tuning done in these models to match observ¬ 
ables. Furthermore, the agreement between the medium 
and high-resolution runs suggests that our results are 
not strongly dependent on numerical resolution, with the 
caveat that we only probe a factor of two in spatial reso¬ 
lution. In the remainder of this paper, we examine how 
the properties of the outflows themselves scale with halo 
mass. 


1 Note that the method used in lAndrews fc Martinil (120131 1 
differs from that of th e other two observed samples. In 
lAndrews &; Martinil (120131 ), the direct method was used to measure 
the metal l icity from stacked spectra of SP SS galaxies, whereas in 
ILee et ahl (120061 ) and ITremonti et al.1 (120041) metallicities were mea¬ 
sured from the flux ratios of strong lines for individual galaxies. 
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Fig. 2.— The cumulative mass loss history for each of the galaxies. In these figures, the mass loss is scaled by the final virial mass of 
the galaxies. The solid black line indicates the total cumulative mass of gas particles accreted onto the disk of the galaxy, including any 
reaccretions of the same particle. The colored solid lines show the cumulative mass of gas particles ever removed from the disk (green), 
ejected such that they become dynamically unbound from the disk (red), and expelled beyond the virial radius (blue). For comparison, the 
dashed line show the total baryonic mass within the disk. Occasionally the dashed line (total disk baryon mass) lies above the black solid 
line (cumulative mass of accreted gas) because of stellar accretion. 
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Fig. 3.— The redshift zero stellar mass-halo mass relation for 
the simulated galaxies (blu e diamonds), compared to abundance 
match ing- derived fits f rom Behroozi et al. (2013), Moste r et al.l 
(2013) and Brook et al. (2014) (grey lines). Filled diamonds repre¬ 
sent high resolution (e = 85 pc) simulations, while empty diamonds 
represent the medium-resolution (e = 170 pc) simulations. In order 
to better mimic the observations, stellar masses were determined 
from mock-photometric observations and halo masses were taken 
from dark matter-only simulations. 



Velocity [km s ’] 


Fig. 4.— The baryonic Tully Fisher Relation for the simu¬ 
lated galaxies (b lue diamonds), compared to observed galaxies from 
IMcGaughl (120051) (grey circles). Filled diamonds represent high res¬ 
olution (e = 85 pc) simulations, while empty diamonds represent 
the medium-resolution (e = 170 pc) simulations. The simulated 
galaxies lie along the same line as the observed galaxies, indicating 
that their baryonic masses scale appropriately with their rotational 
velocities. 

3.2. Baryon fractions in the disk and halo 

The fraction of halo baryons in stars varies strongly 
with halo mass. Dwarf galaxies are known to be con¬ 
siderably less efficient than L* galaxies at forming stars, 
i.e., the mass fraction of the expected halo baryon con¬ 
tent of dwarf galaxies in the form of stars is much smaller 
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Fig. 5.— The z-band Tully-Fisher relation for the simulated 
galaxies (blue diam onds) compared to obs erved galaxies from 
IPizagno et all (120071 4 and IGeha et al.i (2006) (filled grey circles 
and squares, respectively). Filled diamonds show the high reso¬ 
lution sim ulations, empty diam onds the medium-resolution runs. 
While the IPizagno et a □ lf2007f) data are of the t angental rota¬ 
tional velocities, the data from IGeha et all (120061 ) is of HI line 
widths. T herefore, when com paring to the simulated galaxies the 
data from fGeha et all <120060 should be considered lower-limits on 
the actual tangental rotational velocities. In general, the simulated 
galaxies follow the observed Tully-Fisher relation across nearly an 
order of magnitude in velocity. 



Fig. 6 .— The mass-metallicity relation for the simulated galax¬ 
ies (blue diamonds) compared to observed galaxies. Filled dia¬ 
monds represent higher resolution simulations than the empty di- 
amonds. Grey c ircles represent indivi dual observed galaxies from 
I Lee et akl (120061) . Lines show fits from Andrews &; Martinil (2013) 
(dot-dashed line) and ITremonti etTa l. (2004) (the solid line is the 
median of the galaxies in bins of 0.1 dex in stellar mass, the dashed 
line is the contour enclosing 68% of the galaxies, and the dot¬ 
ted li ne is the contour enc losing 95% of the galaxies). All lines 
from [Tremonti et all (120041) were shifted down by 0.26 dex, as the 
method used in that paper produces systematical higher oxygen 
abundances (lErb et al!ll2006l) . 
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Fig. 7. — The top panel shows fraction of the virial mass scaled 
by the cosmic baryon faction in different components of the galaxy. 
In addition to showing the instantaneous galaxy properties at z = 
0 (filled symbols), this plot indicates the total mass ever accreted 
into different components of the galaxy (empty symbols). Blue di¬ 
amonds represent all baryons ever accreted onto the halo (empty 
blue diamonds) and baryons presently within the halo (filled blue 
diamonds). Red squares represent all the baryons ever accreted 
onto the disk (empty red squares) and all baryons presently part of 
the disk (filled red squares) in the form of cold, dense gas (p > 0.1 
amu/cc and T < 1.2 X 10 4 ) or stars. Green stars represent the 
present stellar mass, with values taken direction from the simula¬ 
tion. The bottom panel shows relative efficiency of various methods 
at suppressing star formation. Blue diamonds show the fraction of 
the cosmic baryons accreted to the halo, an indication of the im¬ 
portance of halo preventative feedback, while blue asterisks show 
the fraction of halo gas accreted onto the disk, a result of disk 
preventative feedback. Red triangles show the fraction of baryons 
accreted to the disk that are still there at present day, which is a 
function of the ejective feedback. Green diamonds (mass fraction 
of the disk in the form of stars) show the relative efficiencies of the 
galaxies in forming stars. 

than for L* galaxies. The reason for this must be some 
combination of three factors: preventive feedback, where 
baryons are prevented from accreting either onto the halo 
or onto the disk from the halo; ejective feedback, where 
material enters into the disk but is ejected back into the 
surrounding halo or beyond; and lower global star forma¬ 
tion efficiency, in which gas enters and remains in the 
disk but is less efficient at forming stars. In this section 
we quantify the relative importance of these processes for 
simulated galaxies. 

The top panel of Figure [7] shows the redshift zero bary- 
onic fraction of our simulated galaxies as a function of 
the expected cosmic halo baryon mass (i.e. the baryon 
fraction, /(,, times the halo mass). The filled red squares 
show the z = 0 fraction of baryons in the disk (cold 
gas plus stars), while the filled blue diamonds show the 
fraction of baryons in the halo. The corresponding open 
symbols show the total mass ever accreted into each com¬ 


ponent since a redshift of three (counting a given particle 
only once, even if accreted multiple times). Finally, the 
green stars indicate the present stellar baryon fraction. 

The bottom panel of Figure [7] further summarizes the 
relative importance of preventative feedback, ejective 
feedback, and lower global star formation efficiencies in 
suppressing star formation. In this panel, the fraction of 
the expected baryonic halo content that is accreted onto 
the halo is shown by blue diamonds, and the fraction of 
halo gas accreted onto the disk is shown by blue aster¬ 
isks respectively. The fraction of accreted baryons that 
remain in the disk at z = 0 is shown by red triangles 
and the stellar fraction of the disk is shown by green di¬ 
amonds. As can be seen from the top panel of Figure [71 
all baryon fractions increase with halo mass. The stellar 
baryonic fraction goes from < 1% for Mh a io ;$ 10 1O M 0 
to about 20% for L* halos. 

Preventive feedback manifests in two distinct forms in 
the plot, as shown by the open symbols in the top panel. 
The first is halo preventive feedback, in which the halo 
never receives its cosmic share of baryons. This is seen 
as the difference between the dotted line at unity and 
the open blue symbols. Next is disk preventive feedback, 
whereby the disk does not receive all the baryons ac¬ 
creted onto the halo. This is quantified by the difference 
between the open blue symbols and open red symbols. 

Addressing first halo preventative feedback, at virial 
masses below a few times 10 10 M 0 the halos have re¬ 
duced amounts of material ever accreted onto the halo. 
This halo preventative feedback is also evident in the 
bottom panel where the open blue diamonds mark the 
fraction of the expected baryonic content ever accreted to 
the halo. The existence of halo preventative feedback for 
virial masses less than couple times 10 10 M 0 is expected 
from the impact of the cosmic UV background, which re¬ 
duces baryon f r actions for galaxies with M V j r < 10 9 8 
(|Gnedinl 120001 iHoeft et ~ahl 120061 lOkamoto et al.1 12008D . 
There is some uncertainty in this so-called filtering mass, 
so potentially our simulations are consistent with sim¬ 
ple photoionization suppression, especially since only 
gas that was in the galaxy at z = 3 or later is in¬ 
cluded. However, a distinct effect coul d owe to some 
other form of preventive feedback fe. g.lMo k, Maoll 2002l) 
such as heating by wind ener gy ([Oppenheimer et al.l 
l2010t Ivan de Voort et all l2011al) . At halo masses >~ 
3 x 10 1 o Mq, halos generally accrete their fair share of 
baryons. 

Disk preventive feedback, seen by the difference be¬ 
tween the open blue diamonds and open red squares in 
the upper panel and the blue asterisks in the bottom 
panel, is typically around a factor of two. This has only 
a mild dependence on halo mass, indicating that once 
baryons are inside the halo, an approximately uniform 
fraction of them reach the disk. This trend clearly cannot 
be extrapolated to masses above those considered here, 
since halos well above 10 12 M 0 generally have substantial 
hot gaseous halos but have little gas in their disks (and 
typically do not have disks at all). These simulations are 
also missing possible sources of preventative feedback, 
namely heating from AGN and ionizing radiation from 
the galaxies’ stars, that could possibly in trod uce a mass 
trend within star forming galaxies (jKannan et al.ll2014 ). 
Nevertheless, our simulations here indicate that disk pre- 
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ventive feedback is not strongly dependent on halo mass. 

We now consider ejective feedback, which can be seen 
as the difference between the open and filled points of 
the same color in the top panel. For instance, the halo 
has accreted a baryon fraction indicated by the open blue 
diamonds, but currently only contains the fraction indi¬ 
cated by the filled blue diamonds, so the difference must 
have been ejected from the halo. Similarly, the differ¬ 
ence between the open and filled red squares quantify 
disk ejection. Disk ejection is also apparent from the red 
triangles in the bottom panel, which mark the fraction 
of mass accreted to the disk that remains in the disk at 
z = 0. 

In our simulations, there is only one energetic pro¬ 
cess included that can counteract gravity, eject material, 
and produce galactic winds, namely supernova feedback. 
Hence ejective feedback quantifies how much material 
supernova feedback has removed from the disk and the 
halo. Note that this does not necessarily mean that all 
ejection consists of gas directly heated by the supernova, 
as there could be some entrainment, pushing, or heating 
of the surrounding gas. Ultimately, though, the energy 
source must have been the supernovae. As can be seen 
from the figure, ejective processes have a strong depen¬ 
dence on halo mass. Low-mass disks can eject the vast 
majority of accreted material, whereas for L* galaxies 
the disk ejection is much less. We further explore this 
mass dependency in H3.4I 

Finally, we consider the global star formation effi¬ 
ciency, which is quantified by the difference between the 
solid red squares and the green stars in the top panel 
and by the green diamonds in the bottom panel. Again, 
this shows a significant trend with halo mass, such that 
low-mass galaxies are less efficient at turning their disk 
baryons into stars. In observed galaxies this lower global 
star formation efficiency is evident in the higher gas frac¬ 
tions of dwarf galaxies. The reason is a combination of at 
least two effects. The first is that low-mass galaxies tend 
to be more diffuse, with the lowest mass systems often 
being bulgeless and irregular, and hence the star forma¬ 
tion rate, which is superlinearly dependent on the local 
density, is reduced for a given amount of mass. Second, 
low-mass galaxies have lower metallicity, which reduces 
the amount of dust shielding, and consequentially, the 
H 2 fraction. As a result, the dwarf galaxies have lower 
rates of conversion of disk gas into stars. Overall, ineffi¬ 
ciency in global star formation seem to be of comparable 
strength to ejective effects across all masses. 

So far we have only considered unique accretion events, 
which specifically neglects recycling of previously ejected 
material back into the disk. To examine this we compare 
the disk baryon fractions relative to the total accreted 
and reaccreted fraction onto the disk in Figure [8] In the 
case of the solid red triangles, we consider only unique 
accretions (i.e. not including recycling); this is equivalent 
to red triangles in the bottom panel of Figure [7] In the 
case of the empty red triangles though, we consider all 
accretions, including recycled material. 

Several trends are evident. First, the fraction of both 
unique and total accreted mass that ends up in the disk 
has a strong dependence on halo mass. In the highest 
mass galaxies about 80% of the baryons that are accreted 
once ended up in the disk, while in the lowest masses 
ones it is ~ 4x smaller. The total accretion, in all cases 
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Fig. 8 .— The mass of disk baryons at redshift zero as a fraction 
of the total mass accreted onto the disk. For the filled triangles, the 
denominator is the total mass of unique gas particles accreted. For 
the empty triangles, the denominator includes multiple accretions 
of the same gas particle. According to both definitions, the fraction 
of accreted baryons retained in the disk increases by about a factor 
of four across the halo range covered. 

is also about 4x smaller for unique events than for re¬ 
accreted events. That is, typically there is 3x greater 
recycled accretion than first-time accretion. This number 
is roughly independent of mass, showing that the amount 
of recycling relative to first-time accretion is not strongly 
mas s-de pendent. We explore this topic in greater depth 
in § 13.61 

3.3. Quantifying Ejected Material 

We now examine in more detail the role of ejective 
feedback in regulating the baryonic content for galaxies 
of different masses. Specifically, we quantify the mass of 
outflows that galaxies experience as a function of their 
virial mass. We quantify outflows in several ways: (1) the 
gas that becomes too hot or rarified to be considered part 
of the disk; (2) the gas that becomes energetic enough to 
dynamically escape the baryonic disk (“ejected”); and (3) 
the gas not only ejected from the disk but also expelled 
from the entire halo. Each of these latter categories is a 
subset of the former. We can also subdivide any one of 
those categories between the gas particles which directly 
had supernova energy transferred to them and those that 
were entrained by other affected gas. Note that, with 
regards to Figure [51 all the outflow masses should be 
compared to the total accretion including recycling as 
opposed to unique accretion, since outflows also must 
eject recycled material. 

Figure [9] top panel, shows the total amounts of gas 
ejected from the disk (red squares) and expelled from 
the halo (blue diamonds) since a redshift of 3 as a func¬ 
tion of the virial mass. For comparison, the mass in stars 
formed over the history of the galaxy is also shown (green 
stars). For galaxies with Mh a i 0 < 10 11 M 0 the amount 
of ejected material exceeds the stellar mass and at the 
lowest masses it does so by an order of magnitude. This 
leads to the conclusion that, particularly for smaller ha¬ 
los, the visually-dominant stellar component is merely a 
“trace” amount of mass leftover from the small imbalance 
between accretion and outflows. 

The middle panel of Figure [9] quantifies the mass loss 
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Fig. 9. — In all plots, red squares represent the gas ejected from 
the disk. Top plot: total gas mass ejected over the histories of the 
galaxies as a function of their virial mass. Here the mass of ejected 
gas which eventually crosses the virial radius (blue diamonds) and 
the stellar mass of the galaxies (green stars) are shown for com¬ 
parison. The mass of the outflows is significant compared to the 
stellar mass and increases with galaxy mass. Middle plot: total 
outflow mass divided by the stellar mass formed over the histories 
of the galaxies (the effective mass loading factor) as a function of 
the galaxy’s redshift zero virial mass. In addition to the gas ejected 
from the disk or expelled beyond the virial radius (blue diamonds), 
we also show all the gas that became too hot or rarefied to be con¬ 
sidered part of the disk (filled green circles). The effective mass 
loading factors display strong power law dependencies with mass. 
Bottom plot: mass fraction of gas ever accreted onto the disk of 
the galaxy that is ejected from the galaxy. Empty squares show 
the ejected gas particles that supernova energy was directly trans¬ 
ferred to during the simulation, as opposed to gas particles that 
were entrained. As predicted from the lower baryon fractions of 
dwarf galaxies, higher fractions of gas accreted onto the disk are 
later ejected from less massive galaxies. Lower mass galaxies also 
lose greater fractions of their gas through entrainment rather than 
direct heating. 


more precisely via the effective mass loading factor rj, 
defined as the amount of mass driven out by feedback 
relative to the amount of stellar mass formed. This is 
different than the more canonical mass loading factor 
77 , which is the instantaneous mass loss rate relative to 
the current SFR. The effective mass loading factor is in 
some sense a time-average of 77 over a galaxy’s life. We 
note that, owing to our limited snapshot time resolu¬ 
tion of ~ 100 Myr, fj is actually a lower limit on the 
true value, since there could be recycling happening on 
smaller timescales that would add to the total outflow 
mass. 

The effective mass loading factor, rj , displays the trends 
expected from the top panel. For our most massive halos, 
fj (either for material ejected from the disk or expelled 
beyond the virial radius) is comparable or less than unity, 
showing that galaxies retain in stars roughly as much 
material as they expel. At small masses, however, fj for 
ejected material can reach up to ~ 30, and for the subset 
of the ejecta that is expelled beyond the virial radius, 
fj ~ 10. This relationship demonstrates the extent to 
which less massive galaxies are more efficient at removing 
material from the galaxy through stellar feedback; we will 
quantify this further in the next section. 

We also show the total mass of gas heated by feedback 
relative to the stellar mass as the green circles. Such gas 
is heated enough to no longer be considered part of the 
atomic and molecular ISM but may never have enough 
kinetic energy to be considered “ejected” from the disk. 
Therefore, this group includes gas which is briefly heated 
by supernovae but quickly cools back onto the disk. This 
comprises a significant amount of material and shows the 
impact supernovae can have on the gas available for star 
formation, even when not removing the gas from the disk 
for great lengths of time. 

At large masses, most of the disk baryons are in stars, 
and so fj approximately reflects the amount of ejection 
relative to the disk baryons. But at small masses this 
is not true, because the lower star formation efficiency 
results in large gas reservoirs that far exceed the stellar 
component. Hence a third quantitative view of outflows 
is to compare the outflow mass to the total amount of 
baryons ever accreted onto the disk. 

The bottom panel of Figure [21 shows the fraction of 
ejected material relative to the total accreted mass (in¬ 
cluding recycled accretion). The open squares show the 
ejected material that was directly heated by SNe, while 
the filled squares show the total amount of ejected ma¬ 
terial. As can be anticipated from the reduced baryonic 
fraction of dwarf galaxies, higher fractions of accreted 
material are lost as the virial mass decreases. However, 
the trend is not as strong or as clean as those for fj be¬ 
cause the lower star formation efficiencies of smaller mass 
galaxies results in much of the accretion going into the 
gas reservoir rather than being ejected. 

The effect of entrainment can be seen by comparing 
the open and filled symbols in the bottom panel. At 
high masses, there is little entrainment; most of the gas 
particles either ejected or expelled had energy transferred 
to them directly from supernovae at some point in their 
history. In contrast, there is substantial entrainment of 
gas in dwarf galaxies, in some cases as much as doubling 
the mass of the SN-heated material. 

In order to better analyze the division between all gas 
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Fig. 10.— Fraction of gas heated sufficiently to be removed from 
the disk (Mheated) that was either ejected from the disk (red) or 
expelled beyond the virial radius (blue). The fraction of M heated 
that is ejected strongly decreases with virial mass. This trend 
can be accounted for by the relatively larger disk masses in more 
massive halos. The fraction of M heated that is expelled from the 
halo shows a similar but far weaker trend with mass. Also evident 
is the decreasing fraction of ejected particles later expelled from 
the halo. 

that was heated by supernovae and the subsection that 
was ejected from the disk or expelled beyond the virial 
radius, we show the relative fractions of each for the dif¬ 
ferent halo masses (Figure 1X0)1 . In this figure each gray 
bar represents the total mass of gas that became too hot 
or rarefied to be considered part of the disk (correspond¬ 
ing to the green circles in the middle panel of Figure [9]), 
while the red shows the fraction of that gas that was 
ejected from the disk and the blue the fraction that was 
expelled beyond the virial radius. The fraction of the 
heated gas that was ejected (and, to a lesser extent, ex¬ 
pelled beyond the virial radius) decreases with increasing 
halo mass. This trend demonstrates the greater difficulty 
in gas becoming energetic enough to dynamically leave 
the disk or halo of more massive galaxies. Furthermore, 
the fraction of ejected particles that were later expelled 
beyond the virial radius also decreases with halo mass. 
In the lowest mass galaxies, less than half of the ejected 
gas eventually leaves the halo whereas in the higher mass 
galaxies almost all of the ejected gas does. This trend in 
the fraction of ejected gas later expelled from the halo 
is the result of the relatively larger disks in more mas¬ 
sive halos - the larger mass the disk is relative to the 
halo, the more likely it is that a particle that becomes 
dynamically unbound from the disk will exit the halo. 

3.4. Mass loading factor evolution 

We now examine the mass loading factor, 77 , as a func¬ 
tion of redsliift and halo mass. We obtain a more “instan¬ 
taneous” mass loading factor by taking the total mass 
lost only in the past gigayear at any given redshift and 
comparing it to the mass of stars formed over the same 
time period. While this is not truly instantaneous, it 
allows us to examine how this quantity has varied with 
cosmic time. We also examine 77 as a function of circular 
velocity, as that is more straightforwardly relatable to 
recent models of the physics of outflow driving. In this 
case, “circular velocity” refers to the circular velocity cal¬ 
culated at the virial radius from the total halo mass. 
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Fig. 11.— Total gas mass ejected divided by the stellar mass 
formed within 1 Gyr time bins as a function of halo circular mass. 
The mass loading factors at different redshifts are indicated by 
different color points. A power law fit to all the data points results 
in an exponent of -2.2. 

Figure[lT]shows the mass loading factors for the ejected 
material at four different redshifts, against the corre¬ 
sponding circular velocities at those redshifts. There 
is no discernible trend with redshift in either shape or 
amplitude. Fitting all the values together, we obtain a 
power-law fit of 77 = T7ou c / rc 2 for ejecte d mater ial, roughly 
consi stent with energy-driven winds (jClievalier &; Cleggl 

EMI- 

Cosmological simulations often have to assume a mass 
loading factor since they lack the resolution to directly 
generate outflows. The energy-driven scalings we pre¬ 
dict are consistent with those assumed in state-of-the- 
art simulations that match observed p roperties of galax- 
ies and the circum-galactic medium jDavee1 i _alJj20J_3l 
Ford et, ahl 120131 iVogelsbersfer et aD 12014 iGenel et all 
20141) . They are, however, somewhat different than that 
predicted by the Feedback in Realistic Galaxies (FIRE) 
suite of zoom simulations that also self-consistently drive 
outflows; the FIRE simulations find a shallower depen¬ 
dence for u c i rc > 60 km/s, and a stee per dependence for 
smaller systems (jMuratov et al.llMl5 l). The amplitude is 
somewhat lower than typically assumed in cosmological 
runs; we predict that 77 is unity for u c i rc ~ 100 km/s (i.e. 
Mhaio ~ 10 11 M q ), whereas such simulations typically 
assume unity mass loading at Mh a i 0 ~ 1O 12 M 0 , corre¬ 
sponding to mass loading fa ctors a facto r of two smaller 
than typically assumed (e.g. iDave et al.l(2013 ). However, 
these mass loading factors are only calculated for the 
outflowing gas that is sufficiently energetic to be classi¬ 
fied as “ejected.” If all the gas identified as leaving the 
disk is included, the scaling of the mass loading increases 
to values closer to those assumed in most cosmological 
simulations. 

3.5. Velocities of outflows 

Outflow velocities are gener ally observed to be propor¬ 
tiona l to the circu lar velocity lTMartiiJl2d05t iWeiner et al.1 
120091 though see iSteidel et al.ll2010D . and cosmological 
simulations that employ kinetic feedback typically as¬ 
sume such a scaling. Because our blastwave feedback 
depends on the local density and pressure, affected par- 
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Fig. 12.— Normalized distribution of the velocities of the ejected 
material for individual galaxies. The highest mass galaxies are 
shown in red whereas the lowest mass galaxies are shown in purple. 
The velocities of the ejected material were determined at the step 
immediately following their removal from the disk. The velocities of 
the ejected material are then scaled by the circular velocity defined 
for the halo potential at the time when the particle was ejected. 
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Fig. 13.— Velocities of ejected and expelled material. The me¬ 
dian velocity of the material at the step immediately after ejection 
from the disk is shown as a function of virial mass. These velocities 
are normalized by the circular velocity of the halo at the time of 
ejection. This plot shows both the total ejected gas (red squares) 
and the subsection of that gas which is later expelled (blue di¬ 
amonds). The solid line marks the velocity equal to the circular 
velocity while the dashed line marks the escape velocity of the halo. 

tides show a distribution of energies, temperatures, and 
velocities. We can thus compare these self-consistently 
generated velocities to those observed and employed in 
other models. 

Figure [12] shows the distribution of the velocities of 
ejected gas particles in the snapshot following their re¬ 
moval from the disk scaled by the circular velocity of the 
galaxy at the time of the outflow. Because of the spacing 
between snapshots (100 Myrs), particles are typically al¬ 
ready 0.04 R v i r from the center of the galaxy when their 
velocity is measured. As such, particles may have slowed 
down during the period of time between when they actu¬ 
ally left and when they are identified as having done so. 
Nevertheless, this method allows us to make approximate 
measurements and to compare velocities across galaxies. 
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As in i )3.4l the circular velocity used here is the circular 
velocity calculated at the virial radius for the total mass 
of the halo. Other definitions of rotational velocity will 
produce different scalings. For instance, the asymptotic 
velocities, Vf, used for the baryonic Tully Fisher relation 
(Figured]) are about 1.4 times greater than the circular 
velocities for the halo mass range here. Normalizing the 
outflow velocities by Vf, therefore, would reduce the re¬ 
sults shown by a factor of ~0.7. For the most massive 
galaxies, the median outflow velocities are close to twice 
the circular velocity at the time of the outflow. For the 
least massive galaxies, in contrast, it is clearly smaller. 
This lower velocity of gas ejected from dwarfs galaxies 
is consistent with their relatively smaller fractions of gas 
leaving the halo seen in Figure [TO] Essentially, the rel¬ 
atively smaller disk-to-halo mass ratio of dwarf galaxies 
enables gas to be ejected at relatively smaller velocities. 
However, these smaller velocities are not necessarily suf¬ 
ficient for the ejected material to escape the halo. 

Figure HTT1 quantifies the velocity trend with halo mass 
and circular velocity in more detail. Here we compare the 
median velocities of the ejected material and the subset 
of that material that later escapes the halo at the snap¬ 
shot following their removal from the disk as a function of 
the redshift zero virial mass. Red points show the veloc¬ 
ity of gas ejected from the disk, while blue points show 
the velocity of the subset of that gas that will end up 
expelled from the halo. The ejection velocity increases 
with halo mass and median velocities are between one 
and two and a half times the circular velocity for all halo 
masses. However, there is a huge range of velocities for 
any given galaxy, as indicated by the error bars. Taking a 
closer look, there is a noticeable trend that intermediate- 
mass galaxies have the highest relative outflow velocities. 
There is also a weak trend for material that later escapes 
the halo to have higher velocity, but there is not a strong 
correlation between the velocity at the timestep follow¬ 
ing when it exits the disk and its ability to escape the 
halo. 


3.6. Wind Recycling 

Previous theoretical work has shown that the 
recycling of gas through reaccretion back onto 
the di sk is key to reproducing observed stella r 
masses (iQppenheimer et al.ll20 10; iHenriaues et al.ll20i3h . 
metal enrichment (Daveet al.ll2012D and mass distribu- 
tions dBrook et all 12012‘aj h In this section, we quantify 
re-accretion, or wind recycling, as a function of halo 
mass, as well as the distribution of recycling times. 

In Figure [Tl] we show the fraction of particles that 
has been re-accreted some number of times after being 
ejected from the disk. The number of times particles 
are re-accreted reveals the importance of gas recycling to 
the baryon content of galaxies and the mass of the disks. 
Across all galaxy masses, a significant fraction (20-70%) 
of ejected mass is reaccreted onto the disk of the galaxy. 
There is possibly a slight mass trend towards more mas¬ 
sive galaxies experiencing less reaccretion of their ejected 
material. As in the analysis of the fraction of ejected 
mass that later leaves the halo (Figure HOl). such a trend 
is likely the result of the varying disk-to-halo ratio. Since 
the disk mass is a smaller fraction of the overall halo mass 
for dwarf galaxies, a greater fraction of the particles able 
to escape the disk potential are reaccreted because of the 
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Fig. 14.— Mass fraction of gas ever ejected that was accreted 
some number of times onto the disk as a function of the virial 
mass. Each bar represents an individual galaxy. The red marks 
the fraction of ejected particles never reaccreted during the history 
of the galaxy. The yellow marks the fraction reaccreted once, the 
light green marks the fraction reaccreted twice and so on. Between 
20 and 70% of ejected material is later reaccreted. 

halo potential. 

In contrast, we found no mass trend in the fraction of 
expelled material which is subsequently reaccreted. For 
all galaxies, only about 20% of the expelled material is 
reaccreted at least once. As this process is largely con¬ 
trolled by the growth of the halo (i.e., expelled particles 
can be reaccreted once the potential well is deeper) it 
is reasonable for it to be independent from halo mass 
within a hierarchical growth framework. 

A key parameter discussed in current galaxy formation 
models is the recycling time, i.e. the amount of time 
ejected gas spends outside the disk prior to re-accretion. 
Cosmological simulations predict that the recycling time 
for momentum-driven wind scales roughly inversely with 
the h alo mass for moderate-si zed star forming galax¬ 
ies (iQppenheimer fc Davell2008l l , and some semi-analytic 
models have found that similar scalings are best able to 
reproduc e various observations, such as the stellar mass 
function (IHenriaues et alJfeOlSI) . Here we directly track 
the recycling time in our simulations as a function of halo 
mass. 

Figure [15] shows cumulative histograms of the reac¬ 
cretion time scales normalized by the total number of 
ejected particles. The recycling time distribution has a 
similar shape for all halos. Most gas particles are reac¬ 
creted on short time scales, with most reaccretion taking 
place over 500 Myrs. Since our particle tracking is lim¬ 
ited by our 100 Myr snapshot time resolution, we are 
unable to track particles that have recycling times less 
than 100 Myrs; this fraction is likely not insubstantial. 
Conversely, the recycling time distribution also has a long 
tail, with some gas taking many gigayears to reaccrete. 

In Figure 1151 we also show the median and stan¬ 
dard deviation of the reaccretion times as a function of 
virial mass. This shows a very weak halo mass depen¬ 
dence, roughly f rec ex M^ o , with a typical reaccretion 
timescale of about 1 Gyr. This is in stark contrast to 
cosmological simulations and semi-analytic models that 
seem to favor a strong halo mass depend ence. Inter- 
estingly, the analytic equilibrium model of iMitra et al.1 
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Fig. 15.— Amount of time before reaccretion of ejected particles 
to the disk. Top: normalized cumulative histogram of the amount 
of time particles spend between their ejection and subsequent reac¬ 
cretion onto the disk. The colors represent the virial mass of the 
galaxies, with red being the largest mass halos and purple being 
the smallest mass halos. Where these lines asymptote indicates 
the total fraction of gas ever reaccreted. Bottom: median time for 
reaccretion of ejected particles as a function of halo mass with the 
error bars representing the standard deviation. On average, par¬ 
ticles take 1 Gyr to reaccrete and the recycling times have only a 
very weak dependency with halo mass. 

(l2014f) finds an optimal fit to the stellar mass and halo 
mass, star formation rate, and metallicity from 0 <z <2 
by assuming a recycling time with a weak halo mass de¬ 
pendence of t r ec oc . This is still stronger than the 

dependence our simulations predict, but is closer. 

3.7. Ejecta Properties 

Star formation-driven outflows have been shown to 
have a strong effect on the angular mome ntum distri¬ 
butio n of both baryons and dark matter dBrook et all 
l2012 ai). Outflows preferentially remove low-angular mo¬ 
mentum central gas and, because of cloud-corona interac- 
tions, gas is re-accret ed with higher angular momentum 
(|Marasco et al.ll2012ll . The effectiveness of this processes 
depends on where in the disk most outflow material is 
launched from and where it reaccretes to. Here we ex¬ 
amine these processes over a range of galaxy mass by 
tracking the location and angular momentum of parti¬ 
cles at the time of ejection and reaccretion. 

Figure [16] (top panel) shows the radius within which 
half of the ejected mass originated and is reaccreted as 
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Fig. 16. — Source of the ejected material as a function of virial 
mass. The top panel shows the half-mass radius of the source of the 
ejected material (red squares), the reaccretion of the ejected ma¬ 
terial (green triangles) and the star formation (black stars). Here 
the half-mass radius is defined to be the radius from within which 
50% of the ejected material originates or is reaccreted to or within 
which 50% of the star formation took place. The dotted horizon¬ 
tal lines mark the softening lengths for both resolutions of galaxies. 
The half-mass radii of the source of the ejecta, the location of reac¬ 
cretion, and star formation increase with virial mass, as expected 
from the increasing size of the galaxies. For all galaxies, gas is 
reaccreted to a significantly larger area than it is ejected from. 
The bottom panel normalizes the half mass radius of the ejected 
material by that of the star formation. The source of the ejected 
material roughly follows that of the star formation for galaxies 
above 10 10 Mq. For the four least massive galaxies, the half mass 
radii of star formation is within the softening length (87 pc), which 
limits our ability to draw conclusions as to the relative locations 
of star formation and outflows in these galaxies. Nevertheless, it 
does appear that low mass galaxies eject gas from a broader region 
compared to their star formation. 

a function of halo mass. In our analysis, the originating 
radius is defined by the location of the gas particle at 
the snapshot immediately prior to being ejected and the 
reaccretion radius by the location of the gas particle at 
the snapshot following its reentry into the disk. Similarly, 
the half-mass radius of star formation is the radius within 
which half of the z = 0 stellar mass was formed, which 
we obtain by tracking the location of each star particle 
at the time of creation. 

Both the star formation radius and the ejection radius 
increase strongly with virial mass, roughly following a 
power law R oc which is somewhat shallower than 

the halo virial scaling: R oc M h ' alo . This slight difference 
is consistent with the idea that lower-mass galaxies gen¬ 
erally have a later Hubble type, are more extended, and 
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Fig. 17.— Median angular momentum of the ejected and reac¬ 
creted material as a function of virial mass. The top panel shows 
the median angular momentum of all the ejected material (red 
empty squares) and the ejected material that is later reaccreted 
(red filled squares) at the snapshot immediately prior to it leaving 
the disk. Green triangles indicate the median angular momentum 
of the reaccreted material at the first snapshot after it reenters the 
disk. The bottom panel plots the ratio between the median an¬ 
gular momentum of the material after reaccretion to the median 
angular momentum of it prior to its ejection as a function of virial 
mass. On average, gas is reaccreted at significantly higher angular 
momenta, indicating that it is “spun up” in the halo. 

have more angular momentum. 

In the bottom panel of Figure [TCI we normalized the 
half-mass radius of the ejected material by the half-mass 
radius of star formation. The resulting plot shows that 
the outflowing material roughly follows the star forma¬ 
tion, as expected. There does appear to be a mass trend 
such that more massive galaxies have relatively more con¬ 
centrated outflows. However, this trend is largely set by 
the four lowest mass galaxies, whose star formation half 
mass radii are within the softening length, so it is unclear 
how robust this trend it. It is possible, though, that lower 
levels of star formation Even for the more massive galax¬ 
ies, though, the half-mass radius of ejected material is 
between 0.8 and 2.5 times the SF radius. From the typ¬ 
ically greater half-mass radii of the ejecta, it is apparent 
that while the ejected material is centrally concentrated, 
it is somewhat more dispersed than the star formation. 

Similarly, in Figure 1171 we compare the median an¬ 
gular momentum of the ejected gas immediately prior 
to leaving the disk (open squares), to the median angu¬ 
lar momentum of it immediately after reaccretion (green 
triangles). We also denote the subset of ejected gas that 
is later re-accreted by the filled red squares. There is 
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little distinction between the angular momentum distri¬ 
bution of all the ejected material and the subset of that 
ejected material that is later reaccreted. Hence material 
that recycles does not have a preferential initial angular 
momentum, which implies that it does not come from 
a preferred location in the disk. In contrast, it is clear 
that the reaccreted gas has significantly higher angular 
momentum at the time of reaccretion than it did when 
ejected. 

In the bottom panel of Figure [17] we plot the ratio of 
angular momentum of the gas when reaccreted to that 
when ejected. This figure shows that ejected gas is re¬ 
accreted with higher angular momentum, typically in¬ 
creased by a factor of ~ x2 — 3. There is a large scatter 
and no clear trend with mass, although the highest mass 
galaxies all have fairly low “boosting” of the angular mo¬ 
mentum; a larger sample will be necessary to see if that 
is statistically signific ant. T hese results are consistent 
with previous work bv lBrook et al.l (l2012af) showing that 
ejected gas is “spun up” in the halo before re-accretion. 
It is qualitatively consistent with the redistribution of 
angular momentum necessary to m atch observations o f 
d warf galaxy angular momenta by Ivan den Bosch et all 

(l200l . 

4. DISCUSSION 

The blastwave feedback model used here is essentially 
a variation on the “energy-driven” wind scenario. In 
the energy-driven wind model, gas is driven from galax¬ 
ies assuming SN energy is conserved. This scenario 
is in co ntrast to the “momentum-driven” wind model 
ijMurrav et al.ll2005h in which energy may be dissipated 
away but momentum is conserved. The difference be¬ 
tween the blastwave feedback model and the analytic 
models for energy driven winds applied to many cosmo¬ 
logical simulations is that in the blastwave model the 
transfer of supernova energy to the interstellar media is 
based on the local gas properties and is ignorant of the 
larger galaxy potential. Despite this key difference, the 
scaling for mass loading factor as a function of u c h-c is 
very close to the analytically derived scaling for energy- 
driven winds: u%? c 2 versus v~ irc . We also found that 
the outflow velocity scalings determined for the sample 
of simulated galaxies follow the analytic solution for the 
energy-driven model. The median wind velocities of the 
ejected and expelled material had an approximately lin¬ 
ear scaling with circular velocity. This is consistent with 
the analytic solution for energy-driven winds, namely 
v w = 2<j\JJ l — 1 dMurrav et al.l 1200511 where a is the 
velocity dispersion and, therefore, proportional to the 
circular velocity. 

The correspondence between analytic models and the 
blastwave feedback results imply that the relationship 
between mass loading and galaxy mass comes directly 
from the halo potential. Nevertheless, local gas proper¬ 
ties do affect the efficiency of the energy transfer and, 
there fore, the scaling of the un^ -mass loading relation¬ 
ship. iChristensen et al.l ( 2014b ) showed that the blast- 
wave feedback model results in a much more efficient re¬ 
moval of gas when combined with a highly-resolved inter¬ 
stellar media that included the cold molecular phase, as 
used here. This change in the efficiency is also apparent 
by com paring the results here to those in iWoods et al.1 
(I2014T) . which examined a blastwave feedback model us¬ 


ing lower resolution simu lations and without a molecu¬ 
lar hydrogen model. In IWoods et al.l (I2014T) . the blast- 
wave feedback model was unable to remove baryons from 
ss 10 12 M 0 halos whereas here 10 - 20% of the gas ac¬ 
creted to halos in this mass range were removed by a 
redshift of zero (see Figure [7]). 

A somewhat remarkable aspect of the models shown 
here is that they are able to match the mass-metallicity 
relation despite having energy-driven wind scalings. In 
general, cosmological simulations that include energy- 
driven winds tend to produce a mass-metallicity rela¬ 
tion t hat is too steep (IDave et al.l 1 201. 'll ISomerville et al.1 
1201511 . This may be because such simulations tend to rely 
solely on ejective feedback to suppress star formation at 
low masses (above the filtering mass). The simulations 
presented here, though, show that there is significant pre¬ 
ventive feedback, even above the filtering mass, along 
with a mass-dependent star formation efficiency. For ex¬ 
ample, the fraction of halo gas accreted onto the disk 
is largely independent of halo mass over this sample of 
galaxies (see Figure[7]), indicating that the suppression of 
halo gas accretion onto the disk is similar between dwarf 
and L* galaxies. This preventative feedback can reduce 
the star formation in low-mass galaxies. It does not, how¬ 
ever, have a strong impact on the gas-phase metallicity 
because that is determined by a competition between the 
accretion of fresh gas and the crea tion of new metals from 
stars (e.g. lFinlator fe D ave 2008). Hence it appears that 
to simultaneously match the stellar mass function (or al¬ 
most equivalently, the stellar mass-halo mass relation) 
and the mass-metallicity relation requires having signifi¬ 
cant preventive feedback across all mass scales, not just 
at high masses when active galactic nuclei feedback pu¬ 
tatively happens. In our simulations, this may happen 
owing to wind energy that retards inflow into the disk, 
as also seen bv lvan de Voort et al.l ( 2011b ). 

The energy-driven scaling determined in this suite of 
galaxy simulations follows the general convergence within 
the commun ity towards energy-driven scalings at lowe r 
halo masses (]Ford et al.ll2013t I Vogelsherger et aj~ll20Tll) . 
It should be noted that our model does not include a 
prescription for different forms of early stellar feedback, 
such as radiation pressure, which are generally more con¬ 
sistent with a momentum-driven wind scenario. 

Different models of early stellar feedback have been 
incorporated into a variety of galaxy fo r mation codes 
(ICeverino et al.l 120131: lAumer et all 120131 IStinson et al.1 
120121 iHookins et al.l 1 201. ‘Ib and they have been remark¬ 
able for their ability to reduce early star formation and 
to produce galaxies that match the hig h-redshi f t stel- 
lar mass-halo mass relation (lAgertz fc Kravtsov! 12014 
?). While very similar versions of the code used 
here have been shown to produce dwarf galaxies with 
star formation histories consistent with those observed 
from resol v ed s tellar populations of local dwarf galax¬ 
ies dShen et al.l I2013H . the l arger galaxies likely hav e 
too much early star formation (IChristensen et ah]l2014al f. 
This potential discrepancy suggests that a combination 
of energy-driven and momentum-driven wind may yet 
be necessary to increase outflows sufficiently in higher 
density gas to produce realistic star formation histories. 
An alternative to the addition of early stellar feedback 
exists, however, in improved the modeling o f supernova 
feedback from clusters of young stars, livelier et al.l (|2014f ) 
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describes such a parameter-free, resolution independent 
supernova feedback model. Similarly to early stellar feed¬ 
back models, this model is able to reduce star formation 
before z = 2 and to a produc e Milky Way-mas s galaxy 
with low bulge-to-total mass (iKeller et alJl2015h . 

In addition to affecting the removal of gas, the nu¬ 
merics of the simulation may affect the rates of reac¬ 
cretion of gas. The blastwave feedback model trans¬ 
fers supernova energy directly to the surrounding gas 
particles in the form of thermal energy while simulta¬ 
neously disabling cooling. As a result, gas particles 
may reach very high (> 10 5 K) temperatures while tem¬ 
porarily retaining their high density. The subsequent 
reaccretion time scale will be affected by their ability 
to cool at high temperatures. Other feedback models 
that include less thermal heating of the gas, such as 
those using _a_ h ydrody namicall y-de coupled kinetic wind 
model (ISpringel fe Hernquistll2003afl , should be expected 
to cool and reaccrete more readily. 

Another important numerical factor to consider with 
examining reaccretion rates in simulations is the pres¬ 
ence of thermal instabilities in the halo. As SPH ef¬ 
fectively refines in resolution based on density, the low 
density regions of the halo are resolved to a lesser de¬ 
gree. As a result, instabilities and small scale struc¬ 
ture will be less accurately modeled. Such instabil¬ 
ities become most sensitiv e to num erics when a hot 
gaseous halo is present (lAgertz et al.l 120071 ). However, 
at the masses considered in this study, most simula¬ 
tions sugge st that a virial shock i s not able to be 
supported (iBirnboim &, Dekel 1200. 'll iKeres et~abl I2005L 
though see INelson et al.ll2013l) . 

The version of SPH used in this paper is known to suf¬ 
fer from spurious “surface tension” forces, which make it 
poor at mo deling Kelvin-Helmholtz instabilities. How¬ 
ever, [Hopkins et al. ; 201 20 showed that the changes due 
to different SPH modeling were insignificant compared 
to the feedback models chosen. Furthermore, the net 
amount of gas accretion (as opposed to its thermal state) 
appears to be fairly consi stent across differe nt hydro- 
dynamic implementations (INelson et all [2013 1. Until a 
complete suite of high resolution galaxies can be sim¬ 
ulated with a modified version of the SPH formalism 
(Wadsley et al. 2015 in prep.) or a new hydrodynamics 
methodology ('e.g. lHopkinsll2015l ). these simulations rep¬ 
resent our most complete understanding of the scaling 
outflow properties with mass. 

5. CONCLUSIONS 

Using high-resolution simulations, we have compared 
the outflow properties of twenty galaxy halos spanning 
the mass range from 10 9 5 —10 12 M 0 using particle track¬ 
ing to identify and follow outflowing gas. We show that 
the resulting galaxies match observed global galaxy prop¬ 
erties, indicating that their baryonic content is realistic 
in terms of stellar mass, metallicity, and kinematics. 

Since these trends are determined in a large part by the 
rates of gas inflow and outflow, these simulated galaxies 
present an opportunity to analyze the properties of out¬ 
flows that plausibly produce realistic galaxies. We tracks 
gas particles to identify pristine gas accretion, recycled 
accretion, and gas leaving the disk and/or the halo. From 
this tracking we were able to determine the efficiency of 
outflows, the relative mass of outflowing and recycled 


gas, and the properties of the gas both prior to outflow 
and after accretion. Since the stellar feedback recipe used 
in these simulations depended only on the local proper¬ 
ties of the ISM (not the host galaxy properties), trends in 
the outflow properties with halo mass must result from 
the dynamics of the simulation. 

Our conclusions are summarized as follows: 

• With decreasing galaxy mass, galaxies are signifi¬ 
cantly more efficient at generating outflows. Specif¬ 
ically, mass loading factors show a power-law de¬ 
pendency on circular velocity, with an exponent of 

— 2, which is consistent with energy-driven wind 
models. The similarity between the scaling for the 
analytic solution and that measured for the simula¬ 
tions argues for the mass loading factor being pri¬ 
marily a function of the global halo properties; the 
greater efficiency of small galaxies at driving out¬ 
flows is simply a result of their shallower potential 
well. Furthermore, there was no redshift evolution 
in the scaling of the mass loading factors, which is 
consistent with r\ being primarily a function of halo 
depth. 

• In L* galaxies, ejective feedback is able to reduce 
the baryonic disk mass by 20%, while in galaxies 
with halo masses < 10 10 M 0 there can be as large 
as an 80% mass reduction. These fractions are 
very similar to the stellar-to-disk-mass ratio across 
a range of galaxy masses, indicating that ejective 
feedback is comparable in significance to the glob¬ 
ally averaged star formation efficiency in setting 
the stellar mass of galaxies. Preventative feedback 
also plays an important role in setting the stellar 
mass fraction across the entire mass range; galax¬ 
ies with masses lower than 10 11 M 0 had reduced 
baryon masses compared to the cosmic baryon frac¬ 
tion and across the entire mass range only about 
half of the gas accreted onto the halo was later ac¬ 
creted onto the disk. 

• Recycling is shown to be a common feature of 
galaxy evolution. Approximately 50% of gas 
that becomes dynamically unbound from the disk 
(ejected) is later reaccreted across all masses of 
galaxies. Such recycling occurs primarily on short 
time scales. The median time scales are ~ 1 Gyr, 
with very little dependence on halo mass, and the 
timescales follow a logarithmic fall-off. 

• The source of outflowing material roughly follows 
the spatial distribution of star formation. As 
such, gas is preferentially removed from the cen¬ 
ters of galaxies. Gas is subsequently reaccreted 
with higher angular momenta and farther out in 
the disk, indicating that it is “spun up” through 
interactions with the halo. This trend is consistent 
with previous work indicating that feedback-driven 
outflows can have a significant impact on the an¬ 
gular momentum distribution of disk baryons. 

These results give a quantitative picture for how pre¬ 
ventive feedback, ejective feedback, and star formation 
efficiency plausibly combine to yield the baryon fractions 
in galaxies seen today. In L* halos, the ~ 25% efficiency 
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of baryonic conversion into stars is primarily driven by 
the inability of baryons in the halo to accrete onto the 
disk, presumably owing to t he rapid g row th of a_hot hy¬ 
drostatic halo (e.g. lDekel et, al.ll2009HGabor fc Bournaudl 
120141 ). In modest-sized halos down to ~ 3x 10 10 Mq, ejec- 
tive feedback becomes the dominant modulator of star 
formation, as large-scale cosmological models often as¬ 
sume. At smaller (dwarf) masses, all three effects are 
comparably important: there is significant prevention of 
accretion onto the halo, there is an increasing ejection 
rate, and of the material remaining in the disk, less is 
formed into stars. 

A large number of the outflow properties we exam¬ 
ined are mostly independent of halo mass. Lower mass 
galaxies were more efficient at producing outflows, both 
in terms of their higher mass loading factors and the 
higher fractions of disk gas that was ejected. Indeed, 
the scaling of the mass loading factor with circular ve¬ 
locity calculated from our simulations is consistent with 
previous observations (lArribas et al.1 12014 iRunke et al.1 
1200.4 iChisholm et all [20141) . In contrast, the gas that 
was ejected from the disk tends to have similar v/v c i rc , 
recycling times, and re-accreted fractions across all halo 
masses. These similarities indicate that reducing the 
galaxy mass lowers the threshold for driving outflows, 
but that the outflows themselves are quite similar across 
galaxy mass when scaled by the relevant galaxy prop¬ 
erty. This is somewhat remarkable given the complex 


driving mechanisms and interplay of outflowing gas with 
ambient ISM and halo material. 

Together, our results strongly argue for galactic out¬ 
flows being fundamental to setting the mass of baryons 
and their distribution within galaxies. Metallicities of¬ 
fer yet another strong constraint on the history of the 
baryon cycle. Future work will analyze the source of 
metals, their eventual location, and the relative rates of 
pristine versus enriched accretion. Such metal distribu¬ 
tion mechanisms can be directly constrained by observa¬ 
tions of metal absorption lines in quasar spectra passing 
near gala xies, such as that obtained from the COS -Halos 
program (jTumlinson et al.ll2011k iFord et aLll2015[) . 
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